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Abstract 

The Rouse model has recently been modified to take into account the excluded 
volume interactions that exist between various parts of a polymer chain by incorpo- 
rating a narrow Gaussian repulsive potential between pairs of beads on the Rouse 
chain (Prakash, 2000). The narrow Gaussian potential is characterized by two pa- 
rameters: z* — which accounts for the strength of the interaction, and d* — which 
accounts for the extent of the interaction. In the limit of d* going to zero, the 
narrow Gaussian potential tends to the more commonly used <5-function repulsive 
potential. The influence of the parameter d* , in the limit of infinite chain length, 
on equilibrium and linear viscoelastic properties, and on universal ratios involv- 
ing these properties, is examined here. A renormalization group calculation of the 
end-to-end vector suggests that the value chosen for the variable d* will not affect 
critical exponents, or universal ratios. A similar trend is also observed for results 
obtained with an approximate solution, which is based on the assumption that the 
non-equilibrium configurational distribution function is Gaussian. 
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1 Introduction 



The macroscopic consequences, on both equilibrium and rheological proper- 
ties of dilute polymer solutions, of the microscopic fact that two parts of a 
polymer chain cannot occupy the same place at the same time, have been the 
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subject of innumerable investigations. Analytical treatments of this excluded 
volume effect have, by and large, modeled the repulsive interaction between 
parts of the polymer chain with a Dirac delta potential. As a result, different 
parts of the chain interact with each other only when they are directly in con- 
tact flYamakawa, 1971] ; poi fc Edwards, 1986| ; pttinger, 1989] ; |des Cloizeaux & 



Jannink, 1990| ; |Zylka fc Qttinger, 1991| ; |Schafer, 1999|) . The use of a 5-function 



potential makes it difficult to obtain exact analytical results. Consequently, 
progress has usually been made by using renormalization group methods to 
refine the results of perturbation calculations. On the other hand, numerical 
investigations with Monte Carlo and Brownian dynamics simulations — which 
lead to exact results — have been based on potentials that have a finite range of 
excluded volume interaction, such as the Lennard- Jones potential, because of 
the unsuitability of the ^-function potential for numerical simulations ( praess- 



ley, Hay ward Crest, 1999| ; |Hernandez Cifre Carcia de la Torre, 1999"| ) . The 
quality of the approximate analytical results are usually assessed — at least in 
the case of static properties, which have been much more extensively studied 
than properties away from equilibrium — by comparison with the exact Monte 
Carlo simulations. The comparison is made, in spite of the fact that the an- 
alytical and numerical approaches are based on different potentials, because 
it is generally accepted that the choice of the excluded volume potential does 
not influence the scaling of observables in the limit of long chains. Such com- 
parisons have established that the results of static theories agree well with 
both Monte Carlo simulations and with experimental observations ( |S chafer 
19991 ) . 



Prakash and Ottinger have recently examined the influence of excluded vol- 
ume effects on the rheological properties of dilute polymer solutions by rep- 
resenting the polymer molecule with a Hookean dumbbell model, and using 
a narrow Gaussian repulsive potential to describe the excluded volume inter- 
actions between the beads of the dumbbell (|Prakash fc Ottinger, 199~9| ). The 



same potential has subsequently been used by Prakash ([Prakash, 200(J|) to 



represent the excluded volume interactions between pairs of beads in a bead- 
spring chain model for the polymer molecule. If the bead position vectors 
r u , v — 1, 2, . . . , N, are used to denote the instantaneous configuration of the 
bead-spring chain (which consists of N beads connected together by (N — 1) 
Hookean springs), the narrow Gaussian potential E(r u — r fl ) between any 
pair of beads v and fi can be written as, 



E [r v - Ty) = —3 ) k *T exp 



z* \ . . { H r 2 



UfJ, 



d* 6 J 1 \ 2k B T d** 



where, k-Q is Boltzmann's constant, T is the absolute temperature, H is the 
spring constant, r v/i = r u — r M , is the vector between beads \i and v, and 
the parameters z* and d* are non-dimensional quantities which characterize 
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the narrow Gaussian potential. While z* measures the strength of the ex- 
cluded volume interaction, d* is a measure of the extent of excluded volume 
interaction. Note that in the limit d* tending to zero, the narrow Gaussian po- 
tential becomes a ^-function potential. Compared to the 5-function potential, 
analytical calculations are not significantly harder with the narrow Gaussian 
potential; often, upon setting d* — 0, the predictions of a 5-function poten- 
tial can be obtained. Furthermore, Brownian dynamics simulations can be 
performed with a narrow Gaussian potential. As a result, the use of a nar- 
row Gaussian potential makes it possible to compare approximate analytical 
results and numerical simulations consistently. 

Two approximate solutions have been developed in order to examine the con- 
sequences of using a narrow Gaussian potential (Prakash, 2000| ). The first ap- 
proximation is a first order perturbation expansion in the strength of excluded 
volume interaction, while the second approximation, which is non-perturbative 
in nature, is based on assuming that the configurational distribution function 
is a Gaussian distribution. The equilibrium swelling of the end-to-end vec- 
tor and the radius of gyration, and zero shear rate properties, predicted by 
both the approximate solutions, were compared with the predictions of exact 
Brownian dynamics simulations. 

The Brownian dynamics simulations revealed, rather unexpectedly, that for 
chains with an arbitrary but finite number of beads, the use of a (^-function 
excluded volume potential leads to predictions that are identical to those of 
the Rouse model. In other words theta and good solvents were predicted to 
have the same equilibrium and rheological behavior. Departure from Rouse 
model predictions could only be obtained for non-zero values of d*. As the 
number of beads N became large, however, the simulation results seemed to 
suggest that the use of a 5-function potential may be justified. It must be noted 
though that this conclusion is speculative, and is based on rather limited data. 
It is clearly worthwhile therefore to examine the equilibrium and rheological 
predictions obtained by carrying out the two limits d* — > and N — > oo 
simultaneously. However, it is difficult to obtain data for very long chains 
because of the computational intensity of Brownian dynamics simulations. 

In contrast to Brownian dynamics simulations, both the first order perturba- 
tion expansion and the Gaussian approximation predict a spurious non-trivial 
effect due to the presence of 5-function excluded volume interactions, i.e. when 
d* = 0. In the case of a dumbbell model it was shown rigorously that the source 
of the problem was the incorrect term-by-term integration of a series that was 
not uniformly convergent ( Prakash fc Uttinger, 1999" ). On the other hand, for 



fixed values of N and z*, the Gaussian approximation becomes increasingly ac- 
curate as d* increases. Furthermore, for fixed values of the re-scaled variables 
d = d*/y/~N, and z = z*y/~N, the accuracy of the Gaussian approximation 
increases as N increases ( [Prakash, 2000| ). Since the Gaussian approximation 
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is not as computationally intensive as Brownian dynamics simulations, it is 
possible to obtain predictions in the long chain limit by accumulating data 
for chains of finite length and extrapolating to N — > oo. The Gaussian ap- 
proximation provides a means therefore of exploring the role of the extent of 
excluded volume interaction in the limit of long chains. The purpose of this 
paper is to use the equilibrium and zero shear rate predictions of the Gaus- 
sian approximation to gain insight into the influence of the parameter d* as 
N -> oo. 

As a preliminary motivation to examining the Gaussian approximation in the 
limit of long chains, an elementary renormalization group calculation of the 
equilibrium end-to-end vector is also carried out in this paper. The renor- 
malization group method is a systematic way of examining the macroscopic 
relevance of a microscopic parameter, i.e., if it has macroscopic consequences. 
It is therefore an ideal, albeit technically involved, tool to examine the role of 
the parameter d*. 

The plan of the paper is as follows. The theoretical background for the Gaus- 
sian approximation, and the principal results obtained on making the Gaussian 
approximation are summarized in the next section. The renormalization of the 
equilibrium end-to-end vector is taken up in section 3. The results of extrap- 
olating finite chain Gaussian approximation results to infinite chain length 
are discussed in section 4, and the main conclusions of the present work are 
summarized in section 5. 



2 The Gaussian approximation 



Two fundamental equations, (i) the diffusion equation for the configurational 
distribution function of the polymer chain, and (ii) the Kramers expression 
for the polymer contribution to the stress tensor, are the basic ingredients of a 
kinetic theory of dilute polymer solutions. In the presence of excluded volume 
interactions, one can show that the diffusion equation for a bead-spring chain 
with iV beads, suspended in a Newtonian solvent, has the form flPrakash" 
2000| ), 



dip 
~dt 







N-l 

k B T 7V " 1 
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where, Q { = r i+ i — r i} is the bead connector vector between the beads i and 
i + 1, ip (Qi, • • • , Qjv-i) t) is the configurational distribution function, n(t) is 
the traceless transpose of the velocity-gradient tensor, ( is the bead friction 
coefficient, B kv is an (N — 1) x N matrix defined by, B ku = S k+ i t v — $kv, with 
5 ku denoting the Kronecker delta, and Aj k is the Rouse matrix, 

n [2 for \j - k\ = 0, 

^'fc = Bj V B kv = l—\ f or |j — ^1 = 1, (3) 

I otherwise 



The vector Fl ' is the total excluded volume force on bead v. It is given, in 
terms of the excluded volume potential between the beads of the chain, by the 
expression, 

N Q 

Fi E) = -J2 i-E{r v -r») (4) 

M =l Ul V 



The Kramers expression for the polymer contribution to the stress tensor, for 
a bead-spring chain model with Hookean springs and an arbitrary excluded 
volume potential force, is given by ( [Bird, Curtiss, Armstrong fc Hassager, 
19871 ), 

N-l 

r p = -n p H £ (Q k Q k ) +Z+(N-1) n p k B T 1 (5) 

k=l 



where, angular brackets denote an average performed with the configurational 
distribution function ip (obtained by solving the diffusion equation), n p is the 
number density of polymers, 1 is the unit tensor, and the tensor Z, which 
represents the direct contribution due to excluded volume effects ( Prakasii 
2000Q , is given by, 



N N-l 

II r 

v=l k=l 



p E E B vk {Q k F<P) (6) 



The quantity B vk is a iV x (N — 1) matrix defined by, B vk = k/N — Q(k — u), 
with G (k — v) denoting a Heaviside step function. 

The diffusion equation, Eq. (2), becomes analytically intractable when ex- 
cluded volume interactions are described in terms of the narrow Gaussian 
potential, Eq. (1). As mentioned in the introduction, there are several dif- 
ferent solution schemes, involving either (i) an exact numerical approach, or 
(ii) approximate analytical approaches, that can be developed to overcome 
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this problem. These have been discussed in detail in Prakash (2000|) . Here, we 



summarize one approximate non-perturbative solution procedure, called the 
Gaussian approximation. 

The Gaussian approximation is based on reducing the complex higher order 
moments in Kramers expression, Eq. (5), to functions of only second order 
moments, by assuming that the non-equilibrium configurational distribution 
function ip, is a Gaussian distribution. The Kramers expression can then be 
shown to assume the form, 



N-l 

-p 

k=l 



t- = n»H <*kk + Z + (N - 1) n p k B T 1 (7) 



where, the (N — 1) x (N — 1) matrix of tensor components, cr,- fc = (QjQ k ), is 
the covariance matrix which uniquely characterizes the Gaussian distribution, 
and the tensor Z is now given by the expression, 

1 N 
Z = -z* n p k B T 2J fan ■ n(cr^) (8) 

v, M=l 



Here, the function II (d^) is given by, 

d* 2 l + (T u X l 

n«) = ; L : MJ = (9) 



det ( [d* 2 1 + cr Uf i 



with the tensors ir utl defined by, 

tj max(fi,u) — l 



<V = O^u = <T[iv = 1—^ °jk ( 10 ) 

KbJ- ■ i ■ i \ 



The Gaussian approximation is complete when a scheme for calculating the 
covariance matrix (Tjk is specified. A time evolution equation for can be de- 
rived by multiplying the diffusion equation, Eq. (2), by QjQ k and integrating 
over all configurations. The higher order moments that appear on the right 
hand side of the evolution equation are reduced to second order moments by 
consistently using the ansatz of the Gaussian approximation. The following 
evolution equation for Ojk is then obtained, 

d H JV_1 

~77 Ojk = K ' °3'fe "I" °}'fc ' K ~r \Ojm A m k 4" Aj m CT m k] 

dt S m=l 
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(11) 



where, 




In Eq. (12), the (N — l) x (JV — 1) matrix of tensor components A Jm is defined 
by, 

A?'m = /A {Bj+l,m — B^m) II(<7J + i i/Lt ) — (Bj m — -B^m) Tl{cTjfj) > (13) 

The system of 9 x (JV — l) 2 coupled ordinary differential equations for Oj k , 
Eq. (11), can be solved analytically (correct to first order in velocity gradient), 
by expanding oj fc up to first order in velocity gradient about its isotropic 
equilibrium value. On subsequently making use of Kramers expression, Eq. (7), 
the following first order codeformational memory-integral expansion for the 
polymer contribution to the stress tensor can be derived, 

t 

r p = - J dsG(t-s)j {1] (t,s) (14) 



where, 7m is the codeformational rate-of-strain tensor ( |Bird, Armstrong & 



Hassager, 1987| ), and G(t) is the memory function whose explicit form is given 



in appendix B of Prakash (2000|) . 



Exact expressions for the zero shear rate viscosity, T] Pj o, and the zero shear rate 
first normal stress difference, ^1,0, predicted by the Gaussian approximation, 
can be obtained from Eq. (14). These expressions, which are lengthy, are given 
Prakash (2000|) . It suffices here to note that their evaluation requires the 
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inversion of an (JV — l) 2 x (JV — l) 2 matrix. As a result, the CPU time required 
for their evaluation scales as JV 6 . Generating data for large values of JV becomes 
extremely computationally intensive. Predictions of zero shear rate properties 
have been obtained for chains up to a maximum of JV = 40 beads, since, for 
this value of JV, a single run on an SGI Origin2000 computer with a 195 MHz 
processor required approximately 54 hours of CPU time. 

It is appropriate now to summarize the most significant results obtained on 
making the Gaussian approximation. When the equilibrium swelling and the 
zero shear rate properties, at a constant value of z*, are plotted versus d* for 
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various values of N, and compared with Brownian dynamics simulations, it is 
found that, for each value of N, the Gaussian approximation becomes accu- 
rate beyond a threshold value of d*. However, this threshold value increases 
as N increases, implying that the accuracy of the Gaussian approximation de- 
creases with increasing N. The picture changes considerably, however, when 
the same data is viewed in terms of the re-scaled extent of excluded volume 
interaction d = d*/y/N, and re-scaled strength of the interaction z = z*\fN. 
The Gaussian approximation appears in much better light in the context of 
the re-scaled variables since, at a fixed value of z, it becomes accurate over 
an increasingly larger range of values of d as N increases. Furthermore, as 
clarified below, asymptotic behavior in the limit of large N is observed. 

In the limit of large N, curves for various values of N begin to collapse onto 
a single curve, indicating that all the equilibrium and linear viscoelastic prop- 
erties are independent of the value of N, and are in fact, only functions of 
z and d. This implies that, if some knowledge about the leading order cor- 
rections to the infinite chain length limit can be obtained, data accumulated 
for finite chains for large enough values of N can be efficiently extrapolated 
to the limit N — > oo. As will be seen in the section below, in addition to 
elucidating the role of d* in the theory, a renormalization group calculation of 
the equilibrium end-to-end vector also provides insight into the leading order 
correction. Discussion of the extrapolation procedure, and the results obtained 
on extrapolating the Gaussian approximation data to infinite chains, will be 
taken up in section 4. 



3 Renormalization of the equilibrium end-to-end vector 

The second moment of the end-to-end vector r at equilibrium is given by the 
expression, 



where, the suffix "eq" indicates that the average is carried out with the equi- 
librium distribution function, ip eq , given by, 



The quantity <fi denotes the potential energy of the bead-spring chain, and A/" cq 
denotes the normalization constant. In the present model, the potential energy 
is the sum of the potential energies of all the springs in the chain, and the 




(15) 



j,k=i 




(16) 
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excluded volume interaction energies between all pairs of beads \i and u, 

i N-l i N 

<P=i H T,Qi'Qi+iT, E{r v -r,) (17) 

Z i=l Z M,^ = l 

In order to perform a renormalization group calculation of the end-to-end vec- 
tor, it is first necessary to expand the equilibrium configurational distribution 
function in a perturbation expansion up to first order in the strength of the 
excluded volume interaction. As will be clear shortly, the expansion must be 
made in a space of arbitrary dimensions D. 

Upon expanding both the potential energy (Eq. (17), with E (r v — r M ) given 
by the narrow Gaussian potential, Eq. (1), extended to D dimensions), and 
the normalization constant Af eq , to first order in the strength of the excluded 
volume interaction, one can show that, 



(QiQk) 



eq 



-D 



N 

E 



1 Jdet(d* 2 l + (H/k B T)(r u ,r VfM )f q ) 



°eq, jk 



N 



2 k B T 



Y: (Q 3 Q k E(r u -r,)) 



R 
eq 



where, z* D is an extension of the definition of the z* parameter to D dimensions, 
and angular brackets with superfix U R" and suffix "eq" represent averages 
carried out with the equilibrium distribution function, of the Rouse model, 
extended to D dimensions. As is well known, ip^ is a Gaussian distribution 
function, 



<(Qi,---,Qiv-i; 



H 



27ck B T 



(iV-l)f i 



with a covariance matrix, ^ = (QjQk) eq — {ksT/H) 8j0.. Note that ip eq 
reduces to ip^ in the absence of excluded volume interactions. 

The Gaussian nature of the configurational distribution function in the Rouse 
model has two consequences, (i) the vector r v ^ between beads \i and v is 
also a Gaussian distributed random variable because it is a sum of Gaussian 
variables, and (ii) the complex higher order moment in the second term on 
the right hand side of Eq. (18) can be reduced to second order moments by 
using general decomposition rules for the moments of a Gaussian distribu- 
tion ( [Prakash, 2000| ). On exploiting these consequences, the expression for the 
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covariance matrix (QjQk) eq , correct to first order in z* Dl has the following 
simple form, 



(QiQk) 



eq 



(d* 2 + \fi-i 



(20) 



The function 9 (fx, m, n, v) has been introduced previously in the treatment of 
hydrodynamic interaction ( Pttinger, 1989|) . It is unity if m and n lie between 
fi and v, and zero otherwise, 



9(fi, m, n, v) 



1 if ft < m,n < v or v < m,n < fi 
otherwise 



(21) 



Substituting Eq. (20) into Eq. (15), and carrying out the sum over the indices 
j and k, one can show that the mean square end-to-end vector at equilibrium, 
in a space of arbitrary dimensions D, correct to first order in z* D , is given by, 



(r 2 ) 



cq 



Dk B T 
H 



N 



(N-l) + -z D ^ 

^ (d* 2 + \fi 



(22) 



We now consider the limit of a large number of beads, N. In this limit, the sums 
in Eq. (22) can be replaced by integrals. Introducing the following variables, 



fl v 



(23) 



and exploiting the symmetry in x and y, we obtain, 



(r 2 ) 



eq 



Dk B T 
H 



N{1 + z* D N e/2 dx dy 



[x 





x>y+c 



(d 2 + x — y) 3 e ^ 2 



(24) 



where, e = 4 — D, c is a cutoff parameter of order 1/N which accounts for 
the fact that fi ^ u, and the parameter d = d* / VN, has already been intro- 
duced earlier. It is worth noting that the leading correction to the integrals in 
Eq. (24) is of order N~ l+e l 2 (|Schai'er, 1999 ) — a fact of particular relevance to 
the extrapolation procedure to be adopted in the next section. 

From Eq. (24), it is clear that the excluded volume corrections to the Rouse 
end-to-end vector are proportional to z* D N e l 2 . Therefore, the proper pertur- 
bation parameter to choose is z D = z* D N e l 2 , and not z* D . This well known 
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result of the theory of polymer solutions ( Poi fc Edwards, 1986| |des Cloizeaux 
fc Jannink, 1990) ; ^chafer, 1999 ), indicates that (i) for D = 3 a perturbation 
expansion in z* D is rendered useless for long chains, and (ii) useful results can 
be obtained only when D is close to 4 dimensions. 

The integrals in Eq. (24) can be performed analytically, and the resultant 
expression can be expanded as a power series in e. For the purposes of renor- 
malization, the expansion is required only up to order e°. Carrying out the 
expansion in e, retaining terms up to this order, and neglecting the cutoff c, 
one obtains, 



(r 2 ) 



0(1 



Dk B T 
H 



N 



+ N^ 2 J(d, 



~-K(d,e) 
e 



(25) 



where, 



J(d,e) = (l + 3d 2 ) \(l + d 2 ) e/2 -d 



(26) 



and, 



K{d,e) 



J(d,e) 



[1 + d 2 



-l+e/2 



(27) 



If we keep e finite and set d = in Eq. (25), we get the expression for a 
S-f unction excluded volume potential. This expression has a (1/e) singularity 
which arises because of the neglect of the cutoff. As has been pointed out by 
Ottinger ( pttinger, 1989Q , a renormalization group analysis can be performed 
to get rid of the (1/e) singularity, or, if the cutoff is retained, renormalization 
group analysis can be used to get rid of the cutoff dependence. Curiously, 
when d ^ 0, the parameter d plays a role similar to the cutoff parameter c. 
This follows from the fact that Eq. (25) is free of any singularities as e —>■ 0. 

In order to perform the renormalization of the end-to-end vector, we follow 
standard practice ( pttinger fc Rabin, 1989| ; |Ottinger, 1989| ), and introduce the 
concept of a polymer segment. The polymer segment is used to represent a unit 
of molecular weight, free of any mechanical interpretation. The polymer chain 
is thus assumed to consist of N s segments, each consisting of L s beads, so that 
N = N S L S . The segments are introduced in order to remove the ambiguities 
associated with the arbitrariness in the choice of the number of beads in a 
polymer chain. The size of a polymer segment, L, is assumed to be related 
to the number of beads in a segment L s , through the relation L = Z^L S . 
Introducing a non-dimensional excluded volume parameter, 



(27T 



D_ 

2 z 



D 



U' 2 J(d, e) 



(28) 
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Eq. (25) can be rewritten in terms of N s , u s , and L as, 



<r 2 > 



<-(i 



1 + 



An 2 



(2nN s y /2 Z 



-6/2 



N 



-K(d,e) 



(29) 



The next step in the renormalization procedure consists of introducing a seg- 
ment excluded volume parameter u, which, by renormalizing the bare excluded 
volume parameter u s , takes into account all the excluded volume interactions 
between the many monomers within a segment, 

u = Z u u s (30) 

The parameters u s and L s are then expected to be functions of u, and expanded 
in a Taylor's series, 

u s = u(l + Au + . . .) ; L s = L(l + Bu + . . .) 

so that, for small u, 

Z u = 1 - Au + • • • ; Zjv = 1 - Bu + • • • (31) 

Substituting Eqs. (31) into Eq. (29), keeping only first order terms in u, and 
expanding (27riV s ) e / 2 to first order in e, we have, 

(r% = ^LN. {l + £ g + i In^*) - *M + 2^ B ) } (32) 
Making the choice, 

27r 2 e 4tt 2 

gets rid of all the micro -structure dependent terms in Eq. (32), and leads to,[^] 
^ 2 )e q = ^^Ar s {l + £,ln(2vriV s )} (33) 



The form of Eq. (33) is not consistent with the expectation of a power law 
dependence of (^ 2 ) cq on N s . Renormalization group analysis resolves this prob- 

2 In the limit d — » 0, since K(d, e) is a constant equal to one, the quantity B is 
usually chosen to be equal to (— l/27r 2 e). As a result, an additional term equal to 
(— u/Att 2 ), appears within the braces on the right hand side of Eq. (33). 
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lem by exponentiating the first order perturbation expansion result, giving rise 
to the expression, 



(r\ q = ^L(2n)^N^ (34) 



where, v — ^ + g^2- The exponent v is then obtained by substituting the fixed 
•point value, u = u*, for u ( Pttinger, 1989| ). The fixed point denotes the regime 



where the properties of the infinitely long polymer chain have a power law 
dependence on N s , and the fixed point value is found by setting Z u = 0. For 
instance, in the case of a (^-function excluded volume potential, the fixed point 
value is known to be u* = (7r 2 e/2), which leads to an exponent, v = | + i e. 
In the present instance, the fixed point can be expected to be a function of e 
and d. 

It is not possible to find the fixed point value u* in the present model by 
only considering the renormalization of (r 2 ) eq , since the quantity A does not 
appear in Eq. (32), which is correct only to first order in u. However, we do 
not require the fixed point value in order to make the following argument. The 
form of Eq. (34) suggests that, in the limit N — > oo (which is a consequence of 
Z u = 0), the parameter d*, which always appears as the combination d*/\fN, 
will not in any way alter the N s dependence of the mean square end-to-end 
vector. Indeed, for any finite value of d*, the value of v will always tend, 
as N increases, to the value for a 5-function potential. The validity of the 
above renormalization procedure can only be confirmed by ensuring that the 
same definition of u s and choice of parameter B made here, cancels the micro- 
structure dependence of all the other observable quantities in the theory, such 
as rjp, ^i, etc. This has not been pursued here as it is considered outside the 
scope of the present work. However, the expected independence of the results 
from the choice of d* is in line with the accepted wisdom in static theories of 
polymer solutions that the choice of the excluded volume potential does not 
influence the scaling of observables with molecular weight. Furthermore, di- 
mensionless ratios constructed from observable quantities are also expected to 
be free of micro-structure dependence. In our case, this implies that universal 
ratios can be expected to be independent of the choice of the value of d*. 

We shall make use of the insight gained in this section in our analysis of 
the results of the Gaussian approximation below. This is justified since the 
Gaussian approximation and renormalization group analysis are similar in a 
certain sense; both are exact to first order in the strength of excluded volume 
interaction (the Gaussian approximaton was shown to be exact to first order 
in z* in Prakash (200U| )), and both account for an infinite number of higher 



order contributions. 
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4 The Gaussian approximation in limit of long chains 



C\J 




Fig. 1. Swelling of the radius of gyration versus 1/y/N, at various values of d* . The 
symbols (□ : d* = 0, A : d* = 0.1, : d* = 0.2, O : d* = 0.3, + : d* = 0.4, y : 
d* = 0.5), are the predictions of the Gaussian approximation, while the dashed lines 
are drawn to guide the eye. The filled circle on the y-axis represents the common 
extrapolated value, to N = oo, of all the curves. 



We now examine the aymptotic predictions of the Gaussian approximation, 
in order to examine the role of the parameter d* . The renormalization group 
arguments of the previous section indicate that, for large values of the num- 
ber of beads N, one would expect the scaling with N, of various observable 
quantities, to become independent of d*. A similar behavior is expected of 
non-dimensional ratios constructed with these quantities. 

Figures 1 and 2 examine the dependence of two quantities, predicted by the 
Gaussian approximation, on 1/y/N, for various values of d*, at z = 1. The 
first quantity, a 2 , is both an equilibrium property and a linear viscoelastic 



property, since (as was shown in |Prakash (2000] )), it describes the equilibrium 
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Fig. 2. Non-dimensional ratio constructed from the zero shear rate first normal 
stress difference coefficient and the zero shear rate viscosity (see Eq. (36) for defini- 
tion), versus 1/y/N, at various values of d*. The symbols (□ : d* = 0, A : d* = 0.1, 
: d* = 0.2, Q : d* = 0.3, + : d* = 0.4, \J : d* = 0.5), are the predictions of the 
Gaussian approximation, while the dashed lines are drawn to guide the eye. The 
filled circle on the y-axis represents the common extrapolated value, to N = oo, of 
all the curves. 



swelling of the radius of gyration, and the ratio of the zero shear rate viscosity 
in the presence of excluded volume interactions to the zero shear rate viscosity 
in the Rouse model, 



at. 



(Rl) 



9 'eq 



(35) 



Here, (R 2 g ) is the radius of gyration, and a superfix "i?" on a quantity 
indicates the Rouse model value of the quantity. The second quantity is a 
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Table 1 

Asymptotic values of ratios of equilibrium and zero shear rate properties, at z = 
Numbers in parentheses indicate the uncertainity in the last figure. 



1. 



a 2 a 2 •'' /l " />> 



*i,o/*?o) u *v U R 



1.960 (3) 1.905 (4) 3.32 (1) 0.7324 (4) 0.9723 (3) 



non-dimensional ratio, U<$ v , denned by, 



U Vv = P " lfl (36) 

r lp,0 



The dependence of these quantities on 1/yN was examined for two reasons, 
(i) the leading order correction, in 3 dimensions, to the integrals in Eq. (24) is 
of order 1 / v^A, and (ii) the renormalisation group calculation of the previous 
section, and asymptotic results obtained earlier with the Gaussian approxi- 
mation ( Prakash, 2(J(J(J| )), suggest that the parameter d* always appears in the 



theory as the re-scaled variable d = d* /yN. As a result, the leading order 
corrections to the infinite chain length limit, of all material properties, are ex- 
pected to be functions of 1/yN. In the limit A" — > oo, therefore, all material 
properties should become independent of d*, and depend only on the variable 
z. 

The filled circles on the y-axis in Figs. 1 and 2, represent the common extrap- 
olated value, to A" = oo, of each of the curves for the various values of d*. 
The values for different N, for each value d*, were extrapolated to the limit 
A" — > oo, using a rational function extrapolation algorithm (Press, Tcukolsky. 



Vetterling fc Flannery, 1992| ). Clearly, both a 2 and Uy v become independent 
of d* in the limit A^ — > oo. 

The variables, (i) a 2 = (r 2 ) cq / (r 2 )^, which is the swelling of the end-to-end 
vector at equilibrium, (ii) (^ o/^fo), which is the ratio of the zero shear rate 
first normal stress difference in the presence of excluded volume interactions 
to the zero shear rate first normal stress difference in the Rouse model, and 
(iii) Ur, which is a non-dimensional ratio of the radius of gyration to the 
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end-to-end vector at equilibrium, defined by. 

(r 2 ) eq 



Ur = 6 (37) 



exhibit the same behavior (as that displayed in Figs. 1 and 2), in the limit 
N — > oo. The asymptotic values, for z — 1, of each of these variables is given 
in Table 1.Q 




Z 



Fig. 3. Asymptotic swelling of the radius of gyration versus z. The symbols are the 
results of the Gaussian approximation, while the line is a curve fit using an equation 
of the form given in Eq. (38). The curve fit parameters are given in Table 2. 



Fi gures 3 and 4 describe the asymptotic z dependence of a g and (^1,0/^1 o) 
The behavior of a 2 has not been displayed as it is very similar to that of a 2 



9 'eq 



3 The non-dimensional ratio, U v r = r] p fl/n p r] s (47r/3)(R 

ratio in the present model since it scales with N as N l ~ u . It 
ratio only when hydrodynamic interaction effects are included in the model. 



is not a universal 
aecomes a universal 



17 



120 
100 
o 80 

DC *~ 

* 60 
9? 40 
20 




5 10 15 20 25 30 35 

z 

Fig. 4. Asymptotic ratio of the zero shear rate first normal stress difference coef- 
ficient in the presence of excluded volume interactions to the zero shear rate first 
normal stress difference coefficient in the Rouse model versus z. The symbols are the 
results of the Gaussian approximation, while the line is a curve fit using an equation 
of the form given in Eq. (38). The curve fit parameters are given in Table 2. 

The symbols in Figs. 3 and 4 are the asymptotic predictions of the Gaussian 
approximation. These predictions were obtained, at each value of z, by using 
the extrapolation procedure described above for z = 1. 

The lines through the symbols are curve fits, using an equation which has a 
form commonly used to fit results of renormalization group analysis (|S chafer. 
1999| ), 

y = (l +az + bz 2 )™ (38) 

where, y represents the fitted variable. The parameters a, b and m, for a 2 , 
a 2 and (^o/^fo)' are g iven 111 Table 2. The maximum difference between 
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Table 2 

Parameters appearing in Eq. (38), used to fit the asymptotic predictions of the 
Gaussian approximation, displayed in Figs. 3 and 4. 



a 2 (*i, /*5,) 



a 4.352 4.178 4.058 



b 2.715 2.453 1.755 



m 0.323 0.318 0.630 



the computed data and the curve fit for all the three properties was less than 
0.5%. 

The N dependence of the equilibrium properties (r 2 ) and (R 2 ) , for large 
values of N, can be obtained from Eq. (38), and the definitions of the quantities 
a 2 and a 2 . Both the properties scale identically with molecular weight, namely, 
as N 2u . Experimental results suggest a value of v — 0.592 ±0.003 ( [Hayward &] 
Graessley, 1999| ), while renormalization group calculations and Monte Carlo 
simulations suggest a value of v = 0.588 ( |5chafer, 1999) ; |Graessley et al 



1999] ). Since z scales as iV 1//2 , it is clear from Eq. (38) that for large values 
of N, v — (1 + m)/2. From Table 2, one can see that the values of m for a 2 
and a 2 , imply a Gaussian approximation prediction of v between 0.659 and 
0.662. We expect the value of m for (^ 10 /^f ) to be twice the value of m 
for a 2 . This follows from Eq. (35), since \T/i j0 scales as r] 2 , and r^ and \T/^ 
scale with iV as iV 2 and iV 4 , respectively. As can be seen from Table 2, this 
expectation is reasonably fulfilled. 

It must be noted that the present data, which has been accumulated for rela- 
tively small values of z, might still be describing the crossover region between 
Rouse scaling and the final scaling in the 'excluded volume limit' of large z. In 
that case, the curve fit parameter m might be modified as data is compiled for 
larger values of z. In calculations based on renormalization group arguments, 
the convergence to a value v = 0.588 ± 0.01 is relatively fast (roughly by 
z = 18), but the final convergence is very slow (occuring at z > 100) ( |S chafer j 
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Fig. 5. Universal ratio constructed from the radius of gyration and the end-to-end 
vector versus z. The symbols are the results of the Gaussian approximation, while 
the line is drawn to guide the eye. 



1999). 

The problem with going to larger values of z with the Gaussian approxima- 
tion is that, in order to maintain the accuracy of the numerical extrapolation 
procedure, one must have data for values of N > 40, which, as was men- 
tioned earlier, requires large amounts of CPU time. In the treatment of the 
non-linear microscopic phenomenon of hydrodynamic interaction, the use of 
a normal mode approximation, simultaneously with the Gaussian approxima- 
tion, led to a non-perturbive solution scheme which was as accurate as the 
Gaussian approximation, but significantly less computationally intensive. As 
a result, chains with lengths up to N = 100 could be examined (Prakash & 



Ottinger, 1997| ). It is worth examining if a normal mode approximation in the 
present situation also leads to a significant reduction in the computational 
intensity. 
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Fig. 6. Universal ratio constructed from the zero shear rate first normal stress 
difference coefficient and the zero shear rate viscosity versus z. The symbols are the 
results of the Gaussian approximation, while the line is drawn to guide the eye. 



Figures 5 and 6 examine the asymptotic z dependence of the two universal 
ratios, Ur and U^ v , predicted by the Gaussian approximation. In the Rouse 
model, Ur = 1. Both renormalization group calculations and Monte Carlo 
simulations yield an identical value of Ur = 0.959 in the excluded volume 
limit ( [Schafer, 1999[ ). Though the curve in Fig. 5 decreases relatively rapidly 
for small values of z, and appears to be levelling off as z increases, it still has 
a non-zero slope at z — 30. This suggests that the asymptotic value of the 
ratio, at large values of z, is yet to be reached, and that the present data is 
still in the crossover region. 

For long chains, U^ v = 0.8 in the Rouse model. On the other hand, a renor- 
malization group calculation yields a value U^ v = 0.6288 QOttinger, 1989|) . 
The value of in Fig. 6, at z = 30, is quite close to the prediction of 
the renormalization group calculation (which was based on a ^-function po- 
tential). However, the non-zero slope at z = 30 indicates that the excluded 
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volume limit is yet be reached. 



5 Conclusions 

Renormalisation group arguments, carried out in the context of the end-to- 
end vector at equilibrium, seem to suggest that the scaling with molecular 
weight of the equilibrium and zero shear rate properties, and the values of non- 
dimensional ratios of these quantities, should become independent of d* in the 
limit of large N. The asymptotic results of the Gaussian approximation have 
been examined in the light of these expectations. The data clearly indicates 
independence from the choice of d*. The aymptotic dependence on z presented 
here, of the universal ratios U^^ and Ur, appears to lie in the crossover region 
between Rouse scaling and the scaling in the excluded volume limit of large z. 
Accumulation of more data at larger values of z is required before the scaling 
of the Gaussian approximation in the excluded volume limit can be described 
unambiguously. This task is rendered difficult because of the computational 
intensity of the Gaussian approximation, which scales with chain length as 
A 6 . 
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